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FOREWORD 


Recent  requirements  for  increased  strength  and  service  life  of  machines  and 
structures  have  been  met  by  the  use  of  higher  strength  materials  and  new  fabrication 
and  joining  methods.  Simultaneously,  failures  due  to  fracture  have  increased  rela¬ 
tive  to  those  resulting  from  excessive  deformation.  Frequently  service  conditions 
are  such  that  low  temperature  brittle  fracture,  fatigue  fracture,  and  high  temper¬ 
ature  creep  rupture  must  be  considered  in  a  single  system.  National  concern  with 
increased  safety,  reliability,  and  cost  has  focused  attention  upon  these  problems. 

Methods  are  now  available  to  predict  both  fatigue  crack  initiation  life  and  crack 
propagation  life.  Paradoxically  the  material  properties  required  for  long  fatigue 
crack  initiation  life  are  incompatible  with  the  requirements  of  high  fracture  toughness. 
Thus,  the  conflicting  design  approaches  and  requirements  placed  on  the  material  are 
confusing  and  often  impossible  to  satisfy. 

Numerous  publications  dealing  with  a  variety  of  fracture  problems  have  led  to 
many  new  and  useful  developments.  However,  the  synthesis  of  the  concepts  into 
methods  foi  design,  testing  and  inspection  has  lagged. 

This  program  of  stud/  is  intended  to  contribute  to  the  integration,  correlation, 
and  organization  of  mechanics  and  materials  concepts  and  research  information  into 
a  form  that  will  permit  enlightened  decisions  to  be  made  regarding  fracture  control. 
Reports  are  in  preparation  in  three  categories: 

1.  Research  reports  designed  to  explore,  study  and  integrate  isolated 
and/or  conflicting  concepts  and  methods  dealing  with  life  prediction, 

2.  Reports  to  introduce  and  summarize  the  state-of-the-art  concepts 
and  methods  in  particular  areas,  and 

3.  Example  problems  and  solutions  intended  to  illustrate  the  use  of 
these  concepts  in  decision  making. 

}k  Principal  Investigator 


SUMMARY 


In  a  previous  report,  steady  state  crack  propagation  at  a  constant  stress 
intensity  factor  range  was  analyzed  on  the  basis  of  a  cumulative  fatigue  da  nage 
model.  In  this  model,  the  points  ahead  of  a  crack  tip  were  assumed  to  constitute 
a  series  of  uniaxial  fatigue  specimen  and  fatigue  crack  propagation  was  viewed  as 
the  successive  fatigue  failure  of  these  specimens. 

In  this  report  an  extension  of  the  previous  analysis  to  a  situation  involving  an 
increase  in  the  stress  intensity  factor  range  from  a  lower  level  to  a  higher  level  is 
carried  out.  The  analysis  indicates  that  there  is  a  transient  zone  through  which  the 
crack  has  to  grow  before  attaining  the  full  steady  state  crack  propagation  rate  cor¬ 
responding  to  die  higher  load  level.  Numerical  results  have  been  computed  for  three 
steels  of  comparable  ductility  bui  with  widely  varying  yield  strength.  The  results 
show  that  both  the  crack  length  increment  and  the  number  of  cycles  involved  in 
the  transient  period  decrease  with  increasing  yield  strength.  Also,  for  any  given 
steel,  the  crack  length  increment  during  the  transient  increases  with  the  load  ratio 
but  the  number  of  cycles  is  independent  of  the  load  ratio. 
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INTRODUCTION 


In  a  recent  paper  [1],  steady  state  'atigue  crack  propagation  at  a  constant 
stress  intensity  factor  range  AK  was  analyzed  on  the  basis  of  a  cumulative  fatigue 
damage  model.  In  this  model,  the  points  ahead  of  a  crack  tip  are  assumed  to  con¬ 
stitute  a  series  of  uniaxial  fatigue  specimens,  and  fatigue  crack  propagation  is  viewed 
as  the  successive  fatigue  failure  of  these  specimens.  Such  steady  state  models  have 

also  been  analyzed  by  Fleck  and  Anderson  [2],  Liu  [3],  Rice  [4],  and  McClintock 

[  5J. 

In  this  paper  the  same  model  is  used  to  analyze  the  situation  involving  an 
increase  in  the  stress  intensity  factor  range  from  a  lower  level  to  a  higher  level. 

Such  an  analysis  of  the  transient  behavior  of  fatigue  crack  propagation  appears  to 
not  yet  have  been  carried  out.  Rather,  the  assumption  is  usually  made  that  the 
crack  propagation  rate  after  a  sudden  increase  in  stress  intensity  factor  range 
reaches  the  steady  state  rate  corresponding  to  the  higher  cyclic  load  level  in  only 
a  few  cycles.  The  present  analysis  is  an  attempt  to  verify  this  assumption  and  to 

investigate  the  transient  behavior  as  it  relates  to  the  low  cycle  fatigue  properties 
of  the  material. 

Analysis 

The  problem  to  be  analyzed  is  shown  in  Fig.  la.  A  crack  which  has  been 
propagating  at  a  constant  stress  intensity  factor  range  AK},  is  suddenly  subjected 
to  a  higher  stress  intensity  factor  range  of  AK^  The  X-Y  coordinate  system  is 
stationary  and  has  its  origin  at  the  crack  tip  location  corresponding  to  the  change 
in  loading,  while  the  x-y  coordinate  system  is  attached  to  the  moving  crack  tip.  At 
a  later  instant  the  crack  length  has  increased  by  an  amount  a  (Fig.  lb),  so  that  a 
fatigue  element  at  X  is  at  a  distance  of  x  from  the  current  crack  tip  location  and 
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X  X  +  a.  Initially,  a  -  0  and  X  -  x.  The  initial  damage  of  an  element  at  X  is 
given  by 


00 


D0  <X>  = 


("da/dN)1  J  Nf  (AK^  y) 
X 


dy 


(1) 


where 


da  |  _ 

HnJ  -  steady  state  crack  propagation  rate  at  AK. 


N, 


(AKj,  y)  damage  experienced  in  one  cycle  by  an  element  at  a  distance 


y  from  the  crack  tip  at  stress  intensity  factor  range  AK^. 
The  damage  of  the  same  element  after  N  cycles  is 


N 


D  (X,  N)  =  n0  (X)  +  J  dN 


(2) 


where 


9D 


g^j  ~  damage  per  cycle  = 


Nf(AK2,  x)  Nf  (AK2,  X  -  a) 


Also, 


dN  = 


da 


[da/dN  (a)] 


where 


da 


dN  (a)  "  current  crack  Propagation  rate  when  the  crack  length  increment  is  a. 
If  the  element  at  X  fails  in  N  cycles,  D  (X,  N)  =  1,  and  Eq.  2  reduces  to 


D0  (X)  + 


I 


X 


Nf  (AK2,  X  -  a)  l da/dN  (a)] 


da  =  1 


-■■■' - - - -  - 


Using  Eq.  1,  the  above  can  be  written  as 


where 


X  X 

J  $(a)  L2  (X  -  a)  da  =  L1  dy 


.  da/dN  (a)] 


(da/dN), 


L2  (X  "  a)  Nf  (MZ  X  -  a) 


L,  (y) 


Nf(AKr  y) 


Equation  3  is  a  Volterra  integral  equation  of  the  first  kind  and  has  to  be  solved  for 
the  unknown  function  (a).  The  kernel  is  the  damage  per  cycle  function  and  can  be 
related  to  fne  low  cycle  fatigue  properties  of  the  material  as  follows  [1]: 


A.  (X  +  a.)  M  for  0  £  X  S  R  . 
1  _  (  1  1  P1 

Li  (X)  =  Nf  (AK.,  X)  "  |0  for  X  >  R 


(i  =  1,  2) 

where  R  .  is  the  reversed  plastic  zone  size  under  the  stress  intensity  factor  range 


AK.,  and 

l 


A.  =  2  [ 4 (1  +  n')  7r  E  of'  ef ' ] AK.2^ 


(b  +~cj 


,  usually  varies  between  1  and  2 
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=  (?)  C0Di  +  p*  = 


e  '  AK." 


(no') 

y  ' 


R  .  = - i _ 

P1  [4(1  +  n')  v?  o-  ,21 

y  J 


The  steady  state  crack  propagation  rates  are  given  by 


da  |  _  1  _  Ai“i'|3+1  „  (1+  Rpi|"(3+  * 

dN J.  "3T"/rn  [1  "  “a-2-  J»  i  = 

1  1  ^  ui 


Taking  Laplace  tr  isforms  of  Eq.  3,  using  Eq.  4  and  solving  for  the  Laplace  trans¬ 
form  of  $  (a),  i.e.  $(p), 


*(p) 


=  £^(ai  "  a2}  a!  P)  -  n-j3+  1.  a1  P  +  R  p) 


i'P  +  li  ff?p  +  R  9  p) 


where 


r(-/3  +  1,  a  p)  -  J  e  y  y  ^  dy  is  the  incomplete  Gamma  function, 
ap 

Approximation 

It  is  difficult  to  either  invert  or  obtain  a  useful  asymptotic  solution  for  Eq.  6. 
A  reasrnable  asymptotic  solution  can  be  developed  if  the  assumption  is  made  that 
Rpi  *  00  *  For  usual  values  of  material  properties  the  ratio  Rp./a.  is  0f  the 
order  of  50  to  100.  if  this  assumption  leads  to  a  solution  where  the  crack  propa¬ 
gation  rate  is  within  a  few  percentage  of  the  steady  state  value  when  the  incremental 
crack  length  equals  the  reversed  plastic  zone  size  Rp2f  it  should  give  a  reasonable 
approximation  to  the  true  solution. 
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hquations  5  and  6  can  then  be  simplified  to 


da 

dN 


A.  a.  ~P+  1 
1  i  i 


-  1) 


(7) 


$(p) 


_  ^1  Al  e^°l  “  a2)  a1  P) 

A2  p  -|l(-ft+l,  a9  p) 


(8) 


Though  it  is  extremely  difficult  to  invert  Eq.  8,  it  is  possible  to  obtain  asymptotic 
solutions  of  $(a)  for  small  and  large  values  of  the  crack  length  extension  a  by 
constructing  the  asymptotic  expansions  of  the  right-hand  side  of  Eq.  8  for  large  and 
small  values  of  the  parameter  p  respectively  and  then  applying  the  theorems  of 
Laplace  transform  theory  [6J.  Thus  it  can  be  shown  that 


■5(a)  -  $  [l  -  0  /iL  -  JL_\  +  ft  (ft  +  1)  //  a 
1  \“l  a2j  2  K 


2  .a,2l 


a, 


21 


R2  2 
P  /  a 


2 


and 


^-]  |  ]  +  °  (a3)  as 


(9) 


$(a)  -  3> 


1  +< 


Otr 


ft  -  1 


a. 


ft  -  1 


On 


+  ft  ~  1  /  a  "2 
5-2?^  a  i 


0 


+  (2  -0) 


o'. 


Q1  a2 


ft  -  2 

P2f-21 

a  i 

+ 

l-J 

+  o  (a  ^  Vs  a  - »  to  (10) 


These  two  asymptotic  solutions  can,  in  fact,  be  extended  to  any  number  of  terms 
desired. 


. 
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RESULTS 

Figures  2-4  show  the  variation  of  the  average  crack  propagation  rate  (normalized 

with  respect  to  the  steady  state  crack  propagation  rate  at  AK2)  with  the  crack  length 

increment  a  (as  a  fraction  of  the  reversed  plastic  zone  size  at  AK2)  for  a  low, 

medium,  and  high  yield  steels  of  comparable  ductility  and  for  two  load  ratios 

AK2/AK1  =  2  and  3.  AKj  was  taken  as  20  ksi\/In~.  in  each  case  and  the  material 

properties  used  for  computation  are  shown  in  Table  1.  The  upper  curves  correspond 

to  the  asymptotic  solution  for  large  values  of  crack  length  increment,  a,  taking  two 

and  three  terms  in  Eq.  10,  and  the  lower  curves  correspond  to  the  asymptotic  solution 

for  small  values  of  crack  length  increment,  a,  taking  two  and  three  terms  in  Eq.  9. 

The  crack  propagation  rate  rapidly  approaches  the  steady  state  value  in  each  case, 

reaching  99%  of  the  full  value  for  the  low  yield  steel,  95%  of  the  full  rate  for  the 

medium  yield  steel  and  about  90%  of  the  full  rate  for  the  high  yield  steel  by  the  time 

the  crack  length  increment  equals  the  reversed  plastic  zone  size  R  .  Hence,  the 

pz 

assumption  in  the  previous  section  (R  - >  oo)  is  reasonable. 

The  amount  by  which  the  crack  grows  initially  in  one  cycle  at  the  lower  load 
level  AKj  is  indicated  in  the  lower  part  of  the  figures,  and  it  is  obvious  that  though 
the  crack  propagation  rate  increases  rapidly  with  the  crack  length  increment,  it  takes 
more  than  one  cycle  to  increase  the  average  crack  propagation  rate  significantly. 

Notice  that  considering  only  three  terms  in  the  asymptotic  solution  for  small  crack 
length  increment  is  sufficient  for  the  low  and  medium  yield  steels,  but  not  as  good 
for  the  high  yield  steel. 

Although  theoretically  the  crack  growth  rate  in  every  case  reaches  the  steady 
state  value  exactly  when  the  crack  length  increment  equals  the  reversed  plastic  zone 
size  Rp2,  the  rapidity  with  which  the  crack  growth  rate  approaches  the  steady  state 
value  as  a  function  of  the  crack  length  increment  is  controlled  by  the  parameter  j3 
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(see  Eq.  10).  The  value  of  /3  for  the  three  steels  vary  from  1.  33  for  the  high  yield 
through  1.  51  for  the  medium  yield  steel  to  1.  72  for  the  low  yield  steel.  Conse¬ 
quently,  the  transient  zone  as  a  fraction  of  the  reversed  plastic  zone  size  is  largest 
for  the  high  yield  steel  and  smallest  for  the  low  yield  stei  ..  However,  since  the  re¬ 
versed  plastic  zone  size  for  rhe  high  yield  steel  is  much  smaller  than  the  low  yield 
steel,  the  actual  transient  zone  size  for  the  low  yield  steel  is  larger  than  the  high 
yield  steel. 

In  order  to  determine  how  the  crack  length  increment  varies  with  each  cycle, 

it  is  necessary  to  have  the  solution  for  cr»ck  propagation  rate  for  all  crack  length 

increments.  The  present  asymptotic  solution  provides  this  information  only  for 

small  and  large  values  of  the  crack  length  increment  as  shown  typically  by  the  solid 

lines  in  Fig.  5.  Here  the  asymptotic  solution  for  large  crack  length  increment  has 

been  adjusted  so  that  the  crack  propagation  rate  is  100%  of  the  full  steady  state  value 

when  the  crack  length  increment  equals  the  reversed  plastic  zone  size  R  These 

p2 

two  asymptotic  solutions  are  then  joined  by  the  dotted  lines,  and  the  area  under  the 
complete  curve  from  0  to  a  gives  an  estimate  of  the  number  of  cycles  required  to 
cause  the  crack  length  increment  a.  Figures  6-8  shows  the  estimated  variation  in 
crack  length  increment  with  the  number  of  cycles  for  the  three  steels.  The  dotted 
lines  in  these  figures  correspond  to  crack  length  increments  if  transient  effects  are 
ignored.  Table  2  shows  the  estimated  number  of  cycles  and  crack  length  increment 
necessary  to  achieve  90%,  95%  and  100%  of  the  full  steady  state  crack  propagation 
rate  in  each  case.  Notice  that  the  crack  length  increment  and  number  of  cycles  re¬ 
quired  to  achieve  any  given  fraction  of  the  full  steady  state  value  vary  inversely  as 
the  yield  strengths  of  the  steel.  Also,  for  any  given  steel,  the  number  of  cycles 
required  to  attain  a  given  percentage  ot  the  full  steady  state  value  is  independent  of 
the  load  ratio,  whereas  the  crack  length  increment  increases  with  the  load  ratio. 
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CONCLUSION 

Upon  increasing  the  stress  intensity  factor  range  and  holding  it  steady,  the 
crack  propagation  rate,  according  to  the  present  model,  does  not  immediately  jump 
to  the  '  alue  corresponding  to  the  higher  stress  intensity  factor  range. 

There  is  a  transient  period,  and  the  number  of  cycles  and  crack  length  incre¬ 
ment  necessary  to  achieve  90%  or  95%  of  the  full  steady  state  crack  growth  rate 
decrease  with  increasing  yield  strength  of  the  steel. 

For  any  steel,  the  crack  length  increment  required  to  attain  90%  or  95%  of 

the  fuT  steady  state  crack  propagation  rate  increases  with  the  load  ratio  AK9/AK  , 

2  1 

whereas  the  number  of  cycles  necessary  to  achieve  the  same  crack  growth  rate  is 
independent  of  the  load  ratio. 

If  the  transient  period  is  completely  ignored  and  the  crack  is  assumed  to 
propagate  at  the  steady  state  rate  starting  from  the  first  cycle  after  the  application 
of  the  higher  load,  the  crack  length,  according  to  this  analysis,  is  always  over¬ 
estimated  by  an  amount  of  the  order  of  the  reversed  plastic  zone  size  at  the  higher 
load  level. 
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Crack  Length  Increment  as  Fraction  of  Low  Yield  Steel 

Reversed  Plastic  Zone  Size ,  a/R ^  AK,  =  20ksiv/ln 
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Fig.  2  Crack  Length  Increment  as  Fraction  of  Reversed  Plastic  Zone  Size  at  AK 
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Crack  Length  Increment  as  Fraction  of 
Reversed  Plastic  Zone  Size  at  AK2 ,  a  /R 
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Fig.4  Crack  Length  Increment  as  Fraction  of  Reversed 
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Fig.5  Method  of  Estimating  Crack  Propagation  Rate  for  intermediate 
Values  of  Crack  Length  Increment 
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Fig.6  Estimate  of  Crack  Length  Increment  versus  Cycles 
for  a  Low  Yield  Steel 


